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1 Introduction 


Mutualism is an important biological interaction in nature. It occurs when one 
species provides some benefit in exchange for some beneht, for example, pollina¬ 
tors and flowering plants, the pollinators obtain floral nectar (and in some cases 
pollen) as a food resource while the plant obtains non-trophic reproductive ben- 
ehts through pollen dispersal and seed production. Another instance is ants and 
aphids, in which the ants obtain honeydew food resources excreted by aphids 
while the aphids obtain increased survival by the non-trophic service of ant de¬ 
fense against natural enemies of the aphids. Lots of authors have discussed these 
models [D El 0 [71 nm Ha El El EH E6] . One of the simplest models is the classical 
Lotka-Volterra two-species mutualism model as follows: 

f x(t) = x(t)(ai - bix(t) + ciy(t)), , . 

1 y(t) = yif) («2 - h 2 y{t) + C 2 x{t)). 

Among various types mutualistic model, we should specially mention the fol¬ 
lowing model which was proposed by May [32] in 1976: 

f x{t) = x{t) (n - - Sixit )), 

\ yif) = yit) (r 2 - - e 2 y{t)), 

where x{t),y{t) denote population densities of each species at time t, ri,Ki,bi,ei 
(i=l, 2) are positive constants, ri,r 2 denote the intrinsic growth rate of species 
x(t),y{t) respectively, Ki is the capability of species x(t) being short of y(t), sim¬ 
ilarly K 2 is the capability of species y{t) being short of x{t). For fll.2|h there are 
three trivial equilibrium points 


El — ( 0 , 0 ), E 2 


ri 


'^1 + 


,0), i?3 = (0. 


r2 


^2 + 


and a unique positive interior equilibrium point E* = {x*,y*) satisfying the fol¬ 
lowing equations 


in- - eixit) = 0, 

1 ^2 - - e2y{t) = 0 , 


(1.3) 
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where E* is globally asymptotically stable. 

In addition, population dynamics is inevitably affected by environmental noises, 
May [22] pointed out the fact that due to environmental fluctuation, the birth rates, 
carrying capacity, and other parameters involved in the model system exhibit ran¬ 
dom fluctuation to a greater or lesser extent. Consequently the equilibrium popu¬ 
lation distribution fluctuates randomly around some average values. Therefore lots 
of authors introduced stochastic perturbation into deterministic models to reveal 
the effect of environmental variability on the population dynamics in mathemat¬ 
ical ecology [21 EH ESI EH ESI EH ESI ESI EH EH ESj- Li and Gao et al took into 
account the effect of randomly fluctuating environment in [22] , where they consid¬ 
ered white noise to each equation of the problem fll.2|) . Suppose that parameter 
Vi is stochastically perturbed, with 

ri^ri + aiWi{t), i = 1,2, 

where W\{t), W 2 {t) are mutually independent Brownian motion, z = 1, 2 repre¬ 
sent the intensities of the white noise. Then the corresponding deterministic model 
system fll.2p may be described by the ltd problems: 

f dx{t) = x{t){n - - £ix{t))dt + aix{t)dWi{t), 

\ dy{t) = y{t){r 2 - - £ 2 y{t))dt + a 2 y{t)dW 2 {t). 

On the other hand, population systems may suffer abrupt environmental per¬ 
turbations, such as epidemics, earthquakes, hurricanes, etc. As a consequence, 
these systems are very complex and their sample paths may not be continuous, 
which yields the system fll.dp fail to cope with them. It is recognized that intro¬ 
ducing Levy noise into the underlying population system may be quite suitable 
to describe such discontinuous systems. There exists some interesting literatures 
concerned with SDEs with jumps. We here only mention Bao et al [Hi, Liu 
and Wang [28], Liu and Liang [22]. Motivated by those studies, in this paper we 
consider the following non-autonomous system with jumps: 
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dx{t) = x{t ) [(ri(t) - ~ + ai{t)dWi{t) 

^ + !^li{t,u)N{dt,du)\, , , 

dy(t) = y(t-) [(r 2 (t) - Klit)l%) “ ^ 2 (t)y(t))dt + a 2 {t)dW 2 {t) 

, + j-^^2{t,u)N{dt,du)\, 

where x{t~) and y{t~) are the left limit of x{t) and y{t) respectively, ri{t),hi{t), 
Ki{t), ai(t),i = 1,2 are all positive, continuous and bounded functions on [0, +oo). 
N is a Poisson random measure with compensator N and characteristic measure y 
on a measurable subset Y of (0, +cxo) with /r(Y) < +C)0, N{dt, du) = N{dt, du) — 
li{du)dt, 7 j : Y X n —)■ M is bounded and continuous with respect to /x, and is 
Q3(Y) X 5^t-measurable, i=l, 2. 

In the next section, the global existence and uniqueness of the positive solution 
to problem fll.4p are proved by using comparison theorem for stochastic equations. 
Sections 3 is devoted to stochastic boundedness. Section 4 deals with stochastic 
permanence. Section 5 discusses the persistence in mean and extinction, sufficient 
conditions of persistence in mean and extinction are obtained. 


Throughout this paper, we let (hi, F, {Ff P) be a complete probability space 
with a hltration satisfying the usual conditions. For convenience, we as¬ 

sume that X{t) = {x{t),y{t)) and |X(f)| = ^Jx‘^{t) + y‘^{t). 1 -|- ^i{t,u) > 0,n G 
Y, X = 1, 2, there exists a constant A; > 0 such that 

/ ln{l + u)) V [/n(l -|- 7 j(f, u))]"^fj.{du) < fc, 

Jy 

/3i(t) = O.ba^^t) + [ ['yi{t,u) - ln{l +-fi{t,u))]y{du), i = 1,2, 


Qiit) = / / ln{l +'ji{s,u))N{ds,du), i = l,2, 

Jo Jy 

/ = inf/(A), / = sup/(f). 

t>0 

We end this section by recalling three definitions which we will use in the 
forthcoming sections. 
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Definition 1.1 [1^ If for any 0 < £ < 1, there is a constant (5(e) >0 such that 
the solution X{t) of fll.Sp satisfies 


limsup P{\X(t)\ < (5} > 1 — e, 

t^OO 

for any initial value (xq, yo) > (0,0), then we say the solution X{t) he stochastically 
ultimate boundedness. 

Definition 1.2 If for arbitrary e G (0,1), there are two positive constants 
:= Ci(^) ond C 2 := C 2 (£^) such that 

liminf P{x{t) < Ci} > 1 — £, liminf P{y{t) < Ci} > 1 — £■ 

t^OO t—^OO 

liminf P{x{t) > C 2 } > 1 — liminf P{y{t) > C 2 } > 1 — £■ 

t^OO >-oo 

Then solution of problem mg) is said to be stochastically permanent. 

Definition 1.3 Ifx{t),y{t) satisfy the following condition 

1 1 /■* 
lim - / x(s)ds > 0, lim - / y(s)ds > 0 a.s. 

t^OO t Jq ^ ^ t^OO t Jq ^ ' 

The problem of fll.5p is said to be persistence in mean. 

2 Existence and uniqueness of the positive solu¬ 
tion 

First, we show that there exists a unique local positive solution of fll.hp . 

Lemma 2.1 For the given positive initial value (xo,2/o); there is r > 0 such that 
problem fll.Sp admits a unigue positive local solution X{f) a.s. fort G [0,r). 

Proof: We first set a change of variables : u{t) = lnx(f),u(f) = In^(f), then 
problem fll.Op deduces to 

du{t) = (ri(t) - Pi{t) - - ei(t)e“W)dt + ai{t)dWi{t) 

^ + 7i(M))iV(dt,du), _ x) 

dv{t) = (r 2 (t) - (l 2 {t) - - e 2 (t)e^W)dt + a 2 {t)dW 2 {t) 

+ J.^ln{l +'j 2 {u))N{dt,du) 
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on t > 0 with initial value n(0) = lnxo,n(0) = Inyo- Obviously, the coefficients of 


fl2.ll) satisfy the local Lipschitz condition, then making use of the theorem [21 El] 


about existence and uniqueness for stochastic differential equation there is a unique 
local solution {u{t),v(t)) on t E [0,r), where r is the explosion time. Hence, by 
Ito’s formula, is a unique positive local solution to problem fll.Sp with 

positive initial value. 

Next we need to prove solution is global, that is r = oo. 

Theorem 2.2 For any positive initial value (xo,yo); there exists a unique global 
positive solution {x{t),y{t)) to problem fll.5l) . which satisfies 


X(t) < x{t) < A(f), 9(t) < y{t) < Q(t), t >0, a.s. 


where A(f), A(f), 9{t) and 0(f) are defined as fl2.4l) . fl2.3p . fl2.7l) and fl2.6p . 

Proof: The reference of m was the main source of inspiration for its proof. 
Because of {x{t),y{t)) is positive, from the first equation of fll.Sp . we can define 
the following problem 


J dA{t) = A{t )[[ri{t) - ei{t)A{t))dt + ai{t)dWi{t) + J^'yi{t,u)N{dt,du)], . . 

I A(0) = xo, ^ ’ 


then 


g/o hi h))<^s+/o “1 (s)dW^ (s)+Qi (i) 



is the unique solution of fl2.2p . and it follows from the comparison theorem for 
stochastic equations that 


x{t) < A(f), t G [0,r), a.s. 


(2.3) 


On the other hand. 


m 


e/ohi(«)-/3i(s))'^«+/o cn(s)ciwii(s)+Qi{t) 


Xfj 
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(2.4) 


is the solution to the problem 

I dX{t) = A(t“)[(ri(t) - + ei{t))X(t))dt + ai{t)dWi{t) 

\ + lY'yi{t,u)N{dt,du)], 

( A(0) = Xq, 

then 


x{t) > X{t), t e [0,r), a.s. (2.5) 

Similarly, we can get 

y{t) <e{t), t e[0,T), a.s, (2.6) 


where 


and, 


m 


g/o (>’2 (s) -fe («) )c*S+/o «2 (s)dW2 {s)+Q2(t) 

J_ _|_ g/o'h2(«)-/32(«))rfn+/p a2(«)rfVF2(n)+Q2(s)g2(s)ds 


y(t) > 0(t), t e [0, t), a.s. 


(2.7) 


where 


9{t) 


gfod‘2is)-02is))ds +/p a2{s)dW2{s)+Q2it) 

J_ I „fQ{r2{u)-l32{u))du+f° a2{u)dW2{u)+Q2{s)(r Jg\ i MflWe 

yo Jo ^ ' K 2 {s) ' 


Combining fl2.3p . fl2.5p . fl2.6p with fl2.7p . we obtain 


X{t) < x{t) < A(t), 9{t) < y{t) < Q{t),t > 0, a.s. 

By Lemma 4.2 in [1], we know that A(t), X(t),Q(t),9{t) will not be exploded in 
any hnite time, it follows from the comparison theorem for stochastic equations 
[T5] that {x{t),y{t)) exists globally. □ 


3 Stochastically ultimate boundedness 

In a population dynamical system, the nonexplosion property is often not good 
enough but the property of ultimate boundedness is more desired. Now, let us 
present a theorem about the stochastically ultimate boundedness of fll.5p for any 
positive initial value. 
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Theorem 3.1 Assume that there exists a constant L{q) > 0 such that 


/ \'yi{s,u)\'^fM{du) < L{q), q> I, i = 1,2. 

Jy 

Then for any positive initial value {xo,yo), the solution X{t) of problem fll.Sp is 
stochastically ultimate boundedness. 

Proof: As the reference of |1] we define a Lyapunov function U{x) = x^. By the 


ltd formula: 

rt 

0 
u 


E{e*U{x)) = U{xo) + E e^[U{x) + qx'^ ^dx + ^q{q — l)x'^ ‘^{dx)'^]ds 

= U{^„) + + g[r.(s) - 

= + /y[(l + 7 i(s,ti))'' - 1 - qji(s, u)]/i(du)]U(x)}ds 
< U{xo) + E /g e^{[-ei(s)x + 1 + gri(s) + 

+ /y[(l + 7 i(s,ti))'' - 1 - q^i{s,u)]p{du)\U{x)}ds. 

If g > 1, we can deduce that there exists constant Li{q) > 0 by assumption such 
that 

[/(a:){[(l + qri{s) + - qei{t)x\ 

+ /y[( 1 + 7i(s, u)Y - 1 - g 7 i(s, u)]^^{du)] < Li{q). 

If 0 < g < 1, using (1 + 7 i(s, u)y — 1 — ^ 71 ( 5 , u) < 0, then we have 

I7(x){[(l + qri{s) + - qei{t)x\ 

+ /y[( 1 + 7i(s, u)Y - 1 - g 7 i(s, u)]fi{du)} 

< I7(x)[l + qri{s) — qei{t)x]. 

Therefore, 


Thus, 


E{YU{x)) < E{U{xo)) + Li(g)(e* - 1). 


limsupL^x'^ < Li{q). (3.1) 

t—^OO 


Similarly, we have 


limsupi?y'^ < ^ 2 ( 5 ). 

t^OO 
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(3.2) 








We now combine fl3.1jl . fl3.2p with the formula [x{t)‘^ + y{t)‘^] ^ < 22 + y{ty] 

to yield 

hmsup£^|X|^ < 22[Li(g) + L 2 {q)] < +cxo. 

t—)-CxD 

By the Chebyshev’s inequality ED and the above inequality we can complete the 
proof. □ 


4 Stochastic permanence 


In the study of population models, stochastic permanence is one of the most in¬ 
teresting and important topics. We will discuss this property by using the method 
as in | 2 Sj in this section. 

Theorem 4.1 //min{fi — / 3 i,r 2 — ^ 2 } > 0, then solution of problem fll.511 is 
stochastically permanent. 


Proof: For a positive constant 0 < ?] < 1, we set a function 


Z(x) = -,V(x) = e^^Z^(x). 

X 


Straightforward computation dV{x) by Ito’s formula shows that 

dV (t) = (t) {- Z"^ (t) [ri (t) - 71 (f , u)iJ,{du) - 

Iy { n(i+Yi(t,u))v - - ^] + Z\t){ei{t) + K^j 0 ^)}dt 


- r]aie^^Z^{t)dWi{t) + e^^Z^{t) [(“ l]N{dt,du) 

< r]e^^LZ^{t)dt - r]aie^^Z^{t)dWi{t) + e^^ ’z^{t) - l]N(dt,du), 


Due to 


ln(l + ji(t,u))p.(du) = lim {— - 

17—>-0+ 2 




Iy v{^ + 'li{t,u))^ 

then when f 1 — /3i > 0 , we can choose a sufficiently small 7 to satisfy 


yidu)}, 


n{t)- 


al{t) 


7 i(f, u)y{du) - { ^+ 


l-(l + 7 i(t,n))^ 


Let us choose A > 0 sufficiently small to satisfy 


-<nW-4i 




Y ??(1 + 7i(Lm))'’ 

1 - (l + 7 i(f,M))^ 


jj{du)'^ > 0. 


Iy 7(1 + 7i(Lm))’' 
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Then, there is a positive constant Li satisfying 

, , . 1 — (1 + 7i(t, m))^x , , s 

n(t) - V - X 


{r]-l)al{t) X 


- > -L,. 
V 


dV{t) < T]e^^LZ'^{t)dt-r]aie^^Z^{t)dWi{t) 




1 


( 1+71 (t,u))’' 


— l] N{dt, du), 


where L := Li + + ■^. Integrating and then taking expectations yields 


E[V{t)] = e^^E{Z\x)) < (-)'^ + - 1). 

Xq X 


Therefore, 


limsup E \— 

t^+oo ^x^{t) 


Similarly, when +2 — /32 > 0, we have 


limsup E \— 


< Uli 

- X ■ 


< nL 

- X ■ 


For arbitrary e G (0,1), choosing C 2 (^) = (^)'' and using Chebyshev inequality, 
we yield the following inequalities, 

< 6} = >^}<^ 


P{y(.t) < C 2 } = > 4 } < 

l/nf) C2 


C2 

( 2 ^ 


Hence, 

then. 


limsup F{a;(t) < C 2 } < limsup P{y(t) < C 2 } < 

t —^“|~00 t —^~ 1“00 


liminf P{x{t) > C 2 } > 1 — liminf P{y{t) > C 2 } > 1 — £■ 

i—^~1“00 i—^~1“00 


Combining Chebyshev’s inequality with fl3.ip . fl3.2p . we can prove that for arbitrary 
e G (0,1), there is a positive constant Ci such that 

liminf P{x{t) < Ci} > 1 — liminf P{x{t) < > 1 — e. 

i —^~ 1“00 t — ^~|“00 


This completes the proof. 


□ 


10 















5 Persistence in mean and extinction 


In the description of population dynamics, it is critical to discuss the property of 
persistence in mean and extinction. First, we give a Lemma using the argument 
as in |251E6] with suitable modihcations. 

Lemma 5.1 Suppose that x{t) G x [0,+cx)), i?+). 

(A) If there exist three positive constants T,r] and rjQ such that 

lnx{t) < rjt — rjo f x{s)ds + f ai{s)dWi{s) + Q{t) i = 1 or 2 

Jo Jo 

for all t >T, then 


lim sup 

t —^“|~00 


fo / V 


< — a.s. 
Vo 


(B) If there exist three positive constants T,ri and po such that 

Inxif) > pt — Pq f x{s)ds + [ ai{s)dWi{s) + Q{t) i = 1 or 2 

Jo Jo 

for all t > T, then 

J^x{s)ds p 

limini — -> — a.s. 

i^+oo t TJq 

Proof: (A) Denote Mi(f) = f* ai(s)dWi(s), Qi(t) = f* fy ln(l+'yi(s, u))N(ds, du), 
then Mi{t), Qi{t),i = 1,2 are real valued local martingales vanishing at t = 0. One 
can see that the quadratic variations of Mi{t) and Qiit) are 

{Mi{t),Mi{t)) = [ ai{s)ds < di^t, 


{Qi{t),Qi{t)) = 



(/n(l + 7 i(s, u))Yii{du)ds < kt, i = 1,2 


0 JY 


where (M, M) is Meyer’s angle bracket process, and 


PM{t) = 


d{M, M){s) 


< max{k, di^} 


ds 


< oo. 


ao (1 + ^)^ ^’^'Joil + sy 

By the strong law of large numbers for local martingales [SU], we have 


f ai(s)dWi(s) Qi{t) -10 

hm — -= 0, lim-= 0, a.s. i = 1,2. 

t—^OO t t^OO t 
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Then for arbitrary e > 0, there exists a Ti > 0 such that for t > Ti 


-et < / ai{s)dWi{s) + Qi{t) < et. 


Set g{t) = f* x(s)ds for all t > T 2 , then we have 


dg 


In— <{g + e)t - gog, t>T = max{Ti, T 2 }. 


That is to say : for t > T, integrating this inequality from T to t, 

we can get 

9{t) < 

Therefore 




lim sup 

t —^~1~CXD 


/o xis)ds 


< lim sup 

^ i—^~|~oo Vot 

Using the arbitrariness of e we have the assertion. 

The proof of {B) is similar to (A). The proof is completed. 
Using Lemma 5.1, we have following theorem. 


do 

ln{ 


Tj+e 


) g + 


g + e 
Vo 


□ 


Theorem 5.2 Suppose that fi > /3j, [i = 1,2), X{t) is the positive solution to fll.5p 
with positive initial value {xq^Hq), then the problem fll.Sp is persistent in mean. 


Proof: The method is similar to [22]. We first deduce 


, Inx(f) , Inuft) 

hmsup- - < 0, hmsup- - < 0. a.s. 

t—>CX) t >00 t 


Making use of Ito’s formula to e^lnx, we deduce 

e^lnx-lnxo = [ e''[lnx{s)+ri{s)—l3i{s)-ei{s)x{s)- ]ds+Li{t)+L 2 {t), 

Jo J<i{s) + y{s) 

where Li{t) = e^ai{s)dWi{s), L 2 {t) = J^ln{l + 'yi{s,u))N{ds,du) are 
martingales with the quadratic forms 

{Li{t),Li{t)) = [ e^Wi(s)(is, 

Jo 


12 














{ln{l + ■ji{s,u)))‘^li{du)ds < k / e^^ds. 

Jo 

By the exponential martingale inequality [31], for any positive constants k, 7 , S, we 
can get that 

r{ sup [Li{t) - > de^^lnk} < k-\ 

OKtK'yk 

it follows from the Borel-Cantelli lemma that for almost all u E fl, there is kQ{oj) 
such that for each k > ko{oj), 

Li{t) < 0.5e“'’'^(Lj(f), Li{t)) + Se^'^lnk, 0 <t < 'jk. 

Hence 

eHnx - Inxo < f* e®[/nx(s) + ri(s) - /?i(s) - £(s)x(s) - 

+ /g e^'^af(s)ds + \ke~'^^ e^^ds + 25e'^^lnk 

= J^e^[lnx{s) + ri(s) - fy('yi(s,u) - ln(l + 'Yi(s,u)))/u(du) 

- s(s)x(s) - - |ai(s)[l - - \k[l - e^-^'^Wds 

+ 25e^^lnk 

< Jq e'*[/na;(s) + ri(s) + J^{\^i{s,u)\ + \ln{l +'yi{s,u))\)^{du) 

— |q;i(s)[1 — — ^k[l — e^~'^’‘]]ds + 25e^'^lnk. 

Obviously, for any 0 < s < jk and a: > 0, there is a constant A which is indepen¬ 
dent of k such that 


{L2{t), L2{t)) — / e 


2s 


lnx{s) + ri{s) + 


(|7i(s,w)| -f |/n(l + 7i(s,'u))|)p(dM) 


-0.5a2(s)[l - - 0.5fc[l - < A. 

Then for 0 < t < 7/c, k > fco(<^), we derive 

e^lnx — InxQ < A[e^ ~ 1] + 25e"'^lnk. 


That is 

lnx{t) < e~^lnxo + A[1 — e“*] -1- 2e~^6e*^lnk. 
Letting t —)■ cxo, we have 

, Inx(t) 

hmsup-< 0. 

t—>-00 t 
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Similarly, we get 


t—)-oo t 

On the other hand, applying ltd’s formula to (12.41) we have: 


dlnX{t) = {ri{t) - /3i{t) - {ei{t) + )X{t))dt + ai{t)dWi{t) 

Ki[t) 


That is 


+ / ln{l +'yi{u))N{dt,du). 

Jy 


lnX{t) = lnx{0) + [ (ri{s) - /3i(g) - (£i(s) + ^)^^\ )X{s))ds 

Jo Ki{s) 


For t > T, we have 


+ f Q;i(s)hFi(s)(is + ei(t). 
Jo 


lnX{t) < {ri -/3i + e)t - {ii + ^) / X{s)ds + / a:i(s)lTi(s)cis + 

Ki Jo Jo 

lnX{t) > {ri — (Ji — e)t — {ei -\—/ A(s)(is + / Q;i(s)VFi(s)hs + 

Ki Jo Jo 

Let £ be sufficiently small such that vi — Pi — e > 0, then applying Lemma 5.1 to 
above two inequalities, we get 

Klin - Pi-e) /o X{s)ds f* X(s)ds Ki(ri - Pi + e) 

-3-^- < hm ml — -< hm sup — - < -^-g-. 

bi + eiKi t t bi + SiKi 

Making use of the arbitrariness of £ we get 

fX{s)ds Ki{fi-Pi-e) 

hm sup - > -G-^-. 

t bi + £iKi 


Then 


InA(f) 

hm sup-> 0, a.s. 

t^OO ^ t ~ 


\nx{t) InA(t) 

lim sup- > lim sup-> U, a.s. 

t—>-oo t t—^oo t 
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To sum up, we have 

Similarly, we yield that 

Integrating the hrst equation of fl2.1jl from 0 to t, we yield 

lo dlnx{s) + + - Pi{s))ds 

+ /q ai(s)diyi(s) + Qi{t) - J^ei{s)x{s)ds. 

Because of bi(s)x(s)ds > ds, we obtain 

j bi(s)x(s)ds >^[ - (Inx(f) - Inxo) + /oVi(s) - (3i{s))ds 
+ Jq ai{s)dWi{s) + Qi{t) - /g ei(s)a;(s)ds]. 


, Inx(f) 
hm -= 0. 

t—>oo t 


lim tih!) = 0. 

t^oo t 


which is 


If* - K 

- / (^i(s) + i^i£i(s))a:(s)(is >—^ — (Inx(f) — Inxo) 

t Jo t 


+ (fi—/3i)f+ / Q;i(s)(iiyi(s) + (5i(f)]. 
do 

Since that limt^oo = 0, limt^oo = 0, and limt^oo 


= 0, we get 


px{s)ds h-h 
hm — -> ^ > 0, a.s. 


t—^oo t 


bi + SiKi 


Similarly, we yield 


J^y{s)ds f 2-^2 


hm ° " " " — > > 0, a.s. 


i—>-oo t 


62 + £2-^2 


This completes the proof. 


□ 


Theorem 5.3 Let X{t) he a positive solution of fll.Sp with positive initial value 
X(0), then 

(v4) If vi < $i,r 2 < /32, then x{t),y{t) be extinction. 

(B) If vi > /9i,r2 < ^ 2 , then y{t) is extinction, x{t) is persistent in mean. 

(C) If Ti < /9i,r2 > ^ 2 , then x{t) is extinction, y{f) is persistent in mean. 
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Proof: We first prove Case (^4) of the theorem. Making use of Ito’s formula to 
Inx, X G [0, +cxd) yields 


lnx{t) — lnx{0) < / (ri(s) —/?i(s))(is + / ai(s)(ihhi(s) + 

Jo Jo 

Because of 


/o «i(g)dWi(g) 

t^OO t 

and ri — /9i < 0, we can deduce 


0 , 


lim 

t—yoo 


Qijt) 

t 


0, a.s. 


lim x{t) = 0, a.s. 

t^OO 


Similarly 


lim y{t) = 0, a.s. 


t^OO 


Case {B). Since that r 2 < h, we have limi_j.ool/(t) = 0, a.s. Then 

lnx{t)— lnx{0) < {ri— /3i)t— f iix{s)ds— f hi ds+ f ai(s)(ihhi(s)+ (5i(t), 

Jo Jo Ki Jo 


lnx(t)— lnx{0) > {ri —I3i)t — ei / x{s)ds — bi 

Jo 

Making use of Lemma 5.1, we obtain 


x(s) 


ds-\- / cxi{s)dWi{s)-\-Qi{t). 


Ki{h - fJi) / IoX{s)ds IoX{s)ds ^ Ki{ri - /3i) 


< lim inf 

eiKi + bi 

Hence, we get 


< lim sup 


< 


t—>oo 


Klii + hi 


, a.s. 


i EiKi + bi 

Case {C). Similar to the arguments in Case (H) and {B), it is easy to hnd that: 
x{t) is extinction, y{t) is persistent in mean, if ri < /3i,f2 > ^ 2 - n 
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